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Abstract 

An extraordinary temperature dependence of the resonant Andreev reflection 
via discrete energy level in a normal-metal / quantum-dot / superconductor 
(N-QD-S) system is predicted theoretically by using Green function technique. 



The width of zero bias conductance peak in N-QD-S is about y + T/j 
and does not exhibit thermal broadening, where and F/j are the coupling 
strength between QD and leads. Considering the intra-dot Coulomb inter- 
action, the Coulomb blockade oscillations conducted by Andreev reflection 
differs dramatically from that in N-QD-N. Instead of thermal broadening, fi- 
nite temperature induces more resonant peaks around the oscillation peaks 
of zero temperature. This effect can be applied to determine the coupling 
strength and QD level spacing in N-QD-S. 
PACS numbers: 74.50.+r, 73.40.Gk, 73.20.Dx, 72.15.Nj. 
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Mesoscopic hybrid normal-metal / superconductor (N/S) systems have been investigated 
intensively in the last decade. Aside from their potential applications, these systems pro- 
vide us an opportunity to combine with the two different quantum coherent behaviors, the 
coherence of quasi-particles in mesoscopic system and the coherence Cooper pairs in super- 
conductors. The key mechanism to connect the above two together is the Andreev reflection 
(AR) process at the N/S interface 0,0]. Many new phenomena and effects involving AR have 
been addressed in recent years (for a review see P] and references therein). Among them, 
normal-metal / quantum-dot / superconductor (N-QD-S) is one of the interesting systems, 
in which resonant AR occurs via discrete energy levels of QD Considering the strong 

intra-dot interaction, the problem of Kondo resonance in N-QD-S were also studied by sev- 
eral authors |f^^TT|. In this paper, we shall predict the unique temperature dependence of 



the resonant AR process and investigate the interplay of Coulomb blockade effect and the 
Andreev tunneling in the N-QD-S systems. 

Temperature usually gives us the impression of averaging and smearing. At finite tem- 
perature, the sharp resonant conductance peak usually has a thermal broadening of the 
order of ksT, but we find this is not true for the case of resonant AR conductance. Let us 
begin by considering a N-QD-S system modeled by the following Hamiltonian, 

H = Hl + Hr + H,ot + Ht , (1) 

kcr 

pa p 

Hdot = ^ EqcIc^ , 
Ht = Y^ (thai^c^ + /i.c) + {tnbl^c^ + h.c) , 

ka per 

where and describe the left N lead and the right S lead, Hdot describes QD with one 
spin-degenerate energy level, and Hj- is the tunneling between them. 

In the regime of ksT <^ A, the zero bias conductance is dominated by the AR process, 
and can be derived by using Green function technique (see 0) as, 

f , r df{uj) 



Gnds — -j^ I du 



rl\G[,{u)\' , (2) 



duj 

where f{uj) is the Fermi distribution function, is the coupling strength (defined later), 
is the 12 element of 2 x 2 Nambu matrix C", depicting the conversion from an electron 
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into a hole. C is defined by 



dt e 



-lU)t I 



'^<{c^(t),c{(0)}> <{cT(t),q(0)}>^ 



\^<{cj(t),cj(0)}> <{c;(t),q(0)}> j 
and can be solved from the Dyson equation, 
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and Tl/r = 2ttNl/ji \ tL/R\ , with N^/ji being the density of states in the left /right lead in 
normal state. For comparison we also write down the zero bias conductance for N-QD-N 
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in which single particle tunneling process dominates. In deriving the above formulas, we 
assume that F^/^ is independent of u, which is reasonable in the calculation of zero bias 
conductance |[12|. In the two limits of ksT <^ Tl,Tr and ksT ^ Fi^,F/j, the integrands of 
Eq.(2) and Eq.(8) can be evaluated, and the results are summarized in the following table 



system 


condition 


conductance 


half peak width 


peak shape 


N-QD-S 
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Fig.l shows the curves of Gnds vs Eq and Gndn '^s Eq at different temperatures for the 
symmetric couphng case Tl = Tr = T. One can see the following features: (1) For T ^ 0, 
Gnds has its maximum 4e^//i at £'0 = while Gndn has 2e^//i, and Gnds is much steeper 
than Gndn which is of the Lorentzian shape. (2) With the increase of temperature, both 
Gnds and Gndn are suppressed, but surprisingly, the conductance peak of Gnds does not 
exhibit thermal broadening while Gndn obviously does. Table (9) tells us that the peak 
width of Gnds is always of the order ^yT\ + Tj^ for either <ti Tl^Tr or ksT ^ F^, Vr. 
Feature (1) has been obtained in [Q, while feature (2) is addressed for the first time in this 
work, to our knowledge. 

Taking into account of that AR is a special two-particle process, feature (2) has a simple 
interpretation shown in Fig. 2. For N-QD-N (Fig.2b), the conductance is dominated by the 
process of single-particle tunneling. At finite temperature, the Fermi surface of leads spreads 
out around the chemical potential in the range of ksT, allowing incident electron has energy 
E in this range. The resonant tunneling occurs when Eq = E, resulting in the broadening 
of conductance peak of the order fc^T. For N-QD-S (Fig.2a), on the contrast, the single- 
particle tunneling is forbidden due to superconducting gap, and AR process dominates. In 
the AR process, an incident electron with energy E (with respect to the superconductor 
chemical potential) picks up another electron with energy —E to form Cooper pair and 
enter S, with a hole reflected in the Fermi sea of N. Since both the incident electron and 
the picked electron pass QD through the discrete energy level Eq, the resonant tunneling 
requires not only Eq = E but also Eq = —E. Therefore, the resonant AR process occurs 
only when Eq = 0, i.e., Eq lines up with the chemical potential of superconductor. Because 
Eq has a small broadening and shifting due to coupling with the leads, the condition Eq = 
becomes to \Eq\ < 0(r) with 0(r) being the order of the coupling strength. Thus, the 
resonant AR causes a bottle-neck of the conductance peak width in N-QD-S, leading to a 
broadening about F instead of ksT. 

Next, we investigate a more realistic model, in which QD has multiple single-particle 
energy levels and intra-dot Coulomb interaction, i.e.. 



where i,j are indices of the single-particle energy levels, a, a, (3 are indices of the spin, and U 
is the strength of Coulomb interaction. For N-QD-N with intra-dot Coulomb interaction, it 
is well known that the zero bias conductance vs the gate voltage exhibits Coulomb blockade 




(10) 
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oscillations (CBO) |T^. For N-QD-S in which one of N electrodes is replaced by S, we also 



expect to see CBO except the conducting mechanism is AR. Taking into account of the 
unusual temperature dependence of the resonant AR process, CBO in N-QD-S may have 
significant change from that in N-QD-N. 

We constrain ourselves to discuss the weak coupling and low temperature regime, i.e., 
r -C ksT ^ A. If the coupling between QD and leads is weak enough, the atomic limit 
solution (r 0) should be a good starting point. Because we study the zero bias 
conductance, i.e., Vl — Vr = O"*", QD can be viewed as in equilibrium, and the current 
driven by the small bias voltage servers as a probe to the state of interacting QD. By 
generalizing the exact solution of the retarded Green function of QD in the atomic limit (see 
the Appendix for details) , we propose the following scheme for the calculation of zero bias 
conductance. Suppose QD has / spin-degenerate energy levels, E^, Ei^, E2-1, E21, ■ ■ ■ Ei^, E^, 
with 2^' configurations of occupation represented by F = {N^, Ni^, N21, N21, ■ ■ ■ N^, Ni^) in 
which Ni„ = or 1. In the weak coupling regime, the retarded Green function of QD can 
be derived by an approximated Dyson equation 

G-= g^+g-E^C , (11) 

in which g'' is the atomic limit solution for an isolated interacting QD and "S^ is the self- 
energy caused by the tunneling between QD and leads. S'' is presented in Eq. (6), and g'' 
can be obtained as 

- ^P{F)g^{F) , (12) 



F 



( y. ^ ^ 

Q 1 

\ a;+_B,x{F)+jO+ / 



in which Edot{F) = Y,ia EiaNi^ + f J2ia^j(3 ^iaNjp is the total energy of QD for the config- 
uration F, and Ei^^F) = Ei^ + U Ylijp^iu ^iP renormalized sing-particle levels of QD. 
And the total conductance through QD can be expressed in terms of C" derived from the 
above equations 0. 

We plot Gnds vs Vg and Gndn vs Vg for different temperatures in Fig.3a-3d and Fig.Se, 
where QD contains two spin-degenerate levels, Ei^ = Eii = Ei = Vg, i?2t = E21 = E2 = 
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Vg + AE, with level spacing AE = 0.1, interacting constant U = 1. For <^ AE, we see 
typical CBO pattern in Fig.3a, with nearly equal spaced peaks at = 0.0, 1.0, 2.1, 3.1. 
With the increase of temperature, peak heights are suppressed , roughly proportional to 
(Notice different scaling in Fig.3a-3d). The most striking feature of CBO in N-QD-S 
is, when fc^T ~ AE, temperature induces more resonant peaks around the original CBO 
peaks, instead of thermal broadening. One can see in Fig.3a-3d that the separation among 
peak groups is about U, while the spacing of peaks within one peak group is and each 
resonant peak has a width about F. Notice that all these properties are independent on 
temperature. 

Qualitatively, this unusual pattern can be understood as follows. Due to the intra-dot 
Coulomb interaction, the original two resonances Ei — Ei'^ — -^ij, and E2 — -£'2t — -^21 
exhibits eight sub-resonances, ^1,^2,^1 + U, E2 + U, Ei + 2U, E2 + 2U, Ei + 3U, E2 + 3U 



(see the appendix and Ref. [|T5| for detail). At zero temperature, only four resonances 
El, El + U, E2 + 2U, E2 + 3U are active, and the other four are Coulomb blockaded. But 
for higher temperature of fc^T ~ AE, the Coulomb blockade effect is partially removed. 
Not only Ei but also E2 are active and contribute to the total conductance. For N-QD-N, 
since resonant peak is broadened by finite temperature, the contributions of Ei and E2 
are indistinguishable and combined to form the first CBO peak around = in Fig.3e. 
(This mechanism was first proposed by Y. Meir et al. to explain the anomalous temperature 
dependence of CBO peak heights in N-QD-N |T^.) For N-QD-S, however, resonant AR peak 
has the width of F instead of /c^T. The first group of peaks around = in Fig. 3c consists 
of three distinguishable contributions of AR: AR via Ei, AR via E2, and AR between Ei and 
E2 P]. Other groups of peaks can be understood similarly. Furthermore, one find groups 
with tiny peaks around Vg = 0.5, 1.5, 2.5 in Fig.3d, which can be attributed to the AR 
process between Ei and Ei + U, Ei and E2 + U, E2 and E2 + U, E2 and Ei + U, etc. In short, 
temperature dependence of CBO in N-QD-S differs dramatically from that in N-QD-N. 

Finally, we would like to make three remarks: (1) The extraordinary temperature de- 
pendence of CBO in N-QD-S provides a new approach to determine the coupling strength 
F ~ a/F| + F^ and QD level spacing AE, which are impossible to do in N-QD-N because 
of thermal broadening. (2) The suggested N-QD-S structure is accessible of the up-date 
nano-technology, either by a metal nano-particle connected to superconducting and normal 
metal electrodes [|1^,0, or in a gate controlled N-2DEG-S structure |T8|-pO|]. The key point 
to observe the predicted unusual temperature dependence of AR is to perform resonant AR 
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process via discrete energy levels. (3) Similar extraordinary temperature dependence of zero 
bias conductance have been observed experimentally either in normal-metal / semiconductor 
/ superconductor junctions [^] or in normal-metal / high-Tc superconductor contacts p2 



Perhaps, the discrete electronic state in our model corresponds to the quantum well state 
and surface bound state in their situations. 

In conclusion, we predict an extraordinary temperature dependence of the resonant AR 
process in N-QD-S, in which the conductance peak does not exhibit thermal broadening. 
We also investigate CBO in N-QD-S, and find that finite temperature induces more resonant 
peaks around the original CBO peaks of zero temperature, rather than gives them thermal 
broadening, which can be applied to determine F and AE in N-QD-S. 

This project was supported by NSFC under Grant No. 10074001. T. H. Lin would also 
like to thank the support from the Visiting Scholar Foundation of State Key Laboratory for 
Mesoscopic Physics in Peking University. 
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APPENDIX 



In this appendix, we present the atomic solution in equihbrium, which is generahzed to 
the weak couphng case and hnear response regime in section II. 

In the atomic hmit, QD is nearly isolated from the leads, and the negligible coupling 
to leads determines the equilibrium / non-equilibrium distribution in QD. Suppose QD has 
multiple discrete energy levels, indexed by i = 1,2, •• - L, (here i contains the spin index), 
and QD in the atomic limit can be described by 

L 

Hdot = ^ Eic\ci + UiUj , (Al) 

i=l i<j 

in which the second term is the intra-dot Coulomb interaction, U = e^/2C is the charging 
energy of QD. 

By using the equation of motion, the retarded Green function of QD can be solved 
exactly. For example, for a three-level QD, Hdat — Yli=i EiC^^Ci -\- U{nin2 + n2n^ + ris^i), 
one obtains 

ij\\r _ ((1 -^2)(1 -^3)) ^ ((l-n2)n3) ^ (n2(l-n3)) ^ (nans) 
^^^il^i// - + uj-E,-U + uj-E,-U + uj-E,-2U " ^^^^ 

{uj contains an infinitesimal imaginary part iO^, for all the retarded Green functions). The 
atomic limit solution have clear physics meaning: due to intra-dot interaction, the single 
particle energy level Ei exhibits several sub-resonances, Ei, Ei + U, Ei + 2U ■ ■ ■ , the weights 
of these sub-resonances are determined by the occupation configuration of other energy 
levels. For the three-level QD, Ei has the weight of ((1 — n2)(l — 713)), i.e., the probability of 
both E2 and £'3 are empty; Ei + U has the weights of ((1 — 712)713) and (712(1 — 713)), i.e., the 
probability of one of E2 and £'3 is occupied and the other empty; Ei + 2U has the weight of 
{^2^1.3), i.e., the probability of both E2 and £3 are occupied. 

The calculation of the correlation functions like (77-177,2) should be done sclf-consistcntly 
when QD is in the non-equilibrium distribution. However, for the equilibrium QD in the 
atomic limit, the calculation is straight forward. In this case, the thermal equilibrium QD 
can be depicted by the density operator p = ^e~'^^'*°', and the observable (O) can be 
calculated by Tr{pO). For the three-level QD, the results are fisted as follows 
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(711712^3) 


1 „-/3(Ei+E2+E3+3C/) 


{ni{l - 7Z2)(1 - 713)) 


1 p-pEi 
Z^ 


((1 - ni)n2n3) 


1 p-P(E2+E3+U) 
Z^ 


((1 - 711)712(1 - ^3)) 


1 P-PE2 
Z^ 


{ni{l - 712)713) 


1 p-f5[Ei+E3+U) 
Z^ 


((1 -7li)(l -712)^3) 


Ip-PEs 
Z^ 


(7li7l2(l - ^3)) 


1 p-l5[Ei+E2+U) 
Z^ 


((l-7li)(l-7l2)(l-n3)) 


z^ 



(A3) 



And the weights in Eq.(A2) can be evaluated with the help of this table, for example, 
(712^3) = ((1 - ni)n2n3) + (711712^3)- 

Generally, QD with L energy levels may have 2^ occupation configurations, represented 
by F = (A^i, A^2 • ■ ■ Nl) with AT- = or 1. One has 

P(F)^(n7n,^ = ^e-^^-(^) (A4) 



rrii = 



ii=l 

Ui for Ni = 1 
(1 - m) for Ni^O 



in which Edot{F) = J2i EiNi + U J2i<j ^i^j is the total energy of QD for the configuration 
F. Define = '^i=i{{ci\cl)Y , it is easy to obtain 



(A5) 



F 

L 



in which Ei{F) = Ei + U'^j_^^Nj is the renormalized sing-particle level. 

In short, in the atomic limit, the retarded Green function is the average over different 
occupation configurations weighted by a thermal factor, and each configuration behaves as 
a set of renormalized single particle levels. 
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FIGURE CAPTIONS 



Fig. 1 The zero bias conductance G vs the resonant level £'0 at different temperature for (a) 
N-QD-S and (b) N-QD-N. Parameters are: Vl^Vr^ 0.1, A = 1 for N-QD-S, A = 
for N-QD-N, ksT — 0.001, 0.1, 0.2, 0.3 corresponding to the decrease of peak heights. 

Fig. 2 Schematic picture for understanding the unusual temperature dependence of the con- 
ductance in Fig.l. For (a) N-QD-S, resonant AR requires Eq lines up with the chemical 
potential of the superconducting lead. For (b) N-QD-N, single-particle tunneling oc- 
curs when the resonant level Eq within the range of ksT. 

Fig. 3 Coulomb blockade oscillations of G vs Vg at different temperatures for N-QD-S in 
Fig.3a-3d and N-QD-N in Fig.3e. QD has two spin degenerate levels, E^ — E^ — 
El = Vg and E^st = ^ E2 ^ Vg + AE, with level spacing AE = 0.1, Coulomb 
interacting strength U — 1. Other parameters are: F = Fi = r2 = 0.002, A = 1 for 
N-QD-S, A = for N-QD-N, keT = 0.01, 0.02, 0.05, 0.10 marked in the plots. 
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Fig.2 
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